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We present a local optimal control strategy to produce desired unitary transformations. Unitary 
transformations are central to all quantum computational algorithms. Many realizations of quantum 
computation use a submanifold of states, comprising the quantum register, coupled by an external 
driving field to a collection of additional mediating excited states. Previous attempts to apply control 
theory to induce unitary transformations on the quantum register, while successful, produced pulses 
that drive the population out of the computational register at intermediate times. Leakage of 
population from the register is undesirable since often the states outside the register are prone to 
decay and decoherence, and populating them causes a decrease in the final fidelity. In this work we 
devise a local optimal control method for achieving target unitary transformations on a quantum 
register, while avoiding intermediate leakage out of the computational submanifold. The technique 
exploits a phase locking of the field to the system such as to eliminate the undesirable excitation. 
This method is then applied to produce an 517 (6) Fourier transform on the vibrational levels of the 
ground electronic state of the Na2 molecule. The emerging mechanism uses two photon resonances 
to create a transformation on the quantum register while blocking one photon resonances to excited 
states. 
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PACS numbers: 32.80.Qk, 03.67.Lx, 02.30.Yy 

In recent years there has been growing interest in the 
possibility of realizing quantum computers. Any such 
implementation must consist of a quantum register com- 
prised of a selection of quantum states on which the com- 
putational operations can be performed. A physical real- 
ization of a quantum computer should therefore be able 
to produce these unitary transformations on the register 
by the use of external driving fields pj . 

The general paradigm of quantum computing is to 
break up every computational operation into a sequence 
of simple unitary operations called quantum gates which 
can be considered the basic building blocks of quantum 
computation As the register size grows the number 
of quantum gates required to construct an arbitrary uni- 
tary operation increases rapidly, and therefore the total 
fidelity, which depends on the accumulation of errors at 
each step, decreases drastically It has been shown 
in El IE Op that Optimal Control theory (OCT) methods 
la- lOt. lid. ITU can be utilized in order to calculate a 
field that will directly induce an arbitrary target trans- 
formation. The OCT approach eliminates the need to 
decompose the unitary operator into fundamental oper- 
ations; however the emerging fields are complicated and 
drive the population out of the computational subman- 
ifold at intermediate times. This is undesirable since in 
physical realizations the excited states are prone to de- 
cay and decoherence, and therefore their population at 
intermediate times causes a decrease in the final fidelity. 

In this paper we use a variant of control theory, which 
we refer to as "local control" , to achieve a target unitary 
transformation on a quantum register with no intermedi- 
ate leakage out of the quantum register. The technique 



exploits a phase locking of the field to the system such 
as to eliminate undesirable excitations. We demonstrate 
this approach by obtaining fields to produce an SU(6) 
FT on a quantum register consisting of a submanifold of 
vibrational levels on the ground electronic state of a di- 
atomic molecule 0, without populating the excited 
electronic vibrational states. The emerging mechanism 
uses two photon resonances to transform the quantum 
register while blocking one photon resonances to excited 
states. To some extent the method can be viewed as a 
systematic generalization of the methods of Monroe et al 
[l4| and S0rensen and Maimer to an arbitrary num- 
ber of qubits and arbitrary unitary transformations. The 
more complicated fields presented here can be realized 
using optical pulse shaping techniques 0, 0, . 

The model Hamiltonian consists of a free part Hq and 
an interaction part V — —£iE(t) — (j) E*(t) controlled by 
an external field E(t) through the dipole operator jl — 
Mij K) (i 1 1 with Mij the coupling strengths between the 
ground and excited states \i) and |j) respectively, 



H = H - fXE(t) - ft E*(t). 



(1) 



The system evolution can be described by a time depen- 
dent unitary transformation U(t), the dynamics of which 
is governed by the Schrodinger equation, 



dU{t) 
dt 



-HU(t), 



(2) 



with the initial condition U(0) = I. In order to eliminate 
the free Hamiltonian motion it is common practice to 
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switch to the interaction picture Hamiltonian 

H^H = exp{^H t}Vexp{-^H t} 
n n 

= -fi(t)E(t)-tf(t)E*(t). (3) 

with fi(t) = expij-HotjfLexpi-jrHot}. 

The desired computation is represented by a certain 
unitary transformation O r , which is defined to operate 
solely on the restricted subspace constituting the quan- 
tum register and which must be obtained by the system 
at the final time T. Denoting P r to be a projection onto 
the restricted subspace we define U r {t) = P r U(t)P r to 
be the portion of the evolution operator acting on the 
restricted computational subspace. The goal is therefore 
to obtain a control field E(t) which will maximize the 
overlap, J = |Tr(Oj[/ r (T))| 2 , between the target O r and 
U r (t) at the final time, subject to the constraint that the 
computational manifold population C = Tr (U}(t)U r (t)) 
remains fixed throughout. This constraint enforces the 
elimination of leakage from the computational manifold 
throughout the evolution. 

We introduce, here, a local control method which con- 
tinuously increases the objective J while simultaneously 
holding the constraint C fixed. At each time step, the 
algorithm constructs a field E(t) that fulfills the two 
conditions. The constraint dC/dt = 0, determines the 
direction of the field vector E{t) in the complex plane. 
The sign of the field, however, remains free and is cho- 
sen such as to make the time derivative of J positive, i.e. 
dJ/dt > 0, ensuring an increase in the objective at the 
next time step. 

We now derive the equations determining the direction 
and magnitude of E(t) at each step. Note first that since 
U r = P r UP r we have 

U r = P r ^-HUP r . 

in 

The equation for the constraint is 

C = —Tr(U}U r ) 

= 2Re^Tr(U}P r ^HUP r )X 
= -~Im{Ti(U$P r u(t)UP r )E(t) 

h 

+TiiU t r P r il t (t)UP r )E*(t)} = (4) 

Noticing that the trace operation is invariant under 
transposition of its argument and that the imaginary part 
changes sign under complex conjugation we can take the 
adjoint of the second term and write 

C = -|lm{Tr([/t P r fi(t)UP r )E(t) 

h 

- Tr (P r fl(t ) P r U r ) E{t ) } 

= -2lm{gE(t)} = (5) 



where g = Tr(U^P r p,(t)UP r - P r W ft(t)P r U r ). One can 
now enforce fulfillment of the constraint, eq. JSJ, by 
choosing the electric field to be in the direction of g* 
with (real) proportionality constant a namely 

E(t) = ae* g , (6) 

with e* = g* /\g\ denoting the direction of g* in the com- 
plex plane. We now use the freedom in choosing the sign 
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FIG. 1: The quantities g* and (/ • g) determine respectively 
the direction and magnitude of E(t) in the complex plane. 

of a to assure monotonic increase in the objective. 

j = ±\Tr(O r U r )\ 2 

= 2Re|Tr(Ot C 7 r )*l T r(OtC7 r )| 

= 2Re |r?*Tr(OjP r i#;7P r )| 

= a2Re{/e;} > 0, (7) 

where we define r\ = Tr{0\Ur) and / = 
--i-Tr {r]*O^P r fi{t)UP r - rjP r W ft(t)P r O r } ■ Note that 
in advancing to the last line we have employed the same 
sequence of steps as in eq. 0} and |[SJ| . In order to enforce 
an increase in J, we choose 

a = Re{fg*} 

= (f'9), (8) 

with (/ • g) denoting a scalar product in the complex 
plane, thus guaranteeing that dJ/dt > 0. The electric 
field is therefore chosen as 

E(t) = (f ■ g)e* g , (9) 

It is important that E(t) depend on the magnitude of 
g* and not only on its direction e*, since for a vanishing 
\g\ the direction e* becomes undefined and numerically 
unstable. We therefore wish that E be proportional to 
\g\ such that a phase jump in g be accompanied by a 
vanishing of E, thus avoiding abrupt phase jumps in the 
emerging field. 
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Note that there is still freedom in determining the mag- 
nitude of E(t), which can be utilized to control the order 
of magnitude of the field strength by multiplying by a 
positive envelope function, S{t). It is sometimes neces- 
sary for numerical reasons to limit the maximum allowed 
field such that \E{t)\ < E max . This can be achieved by 
transforming a E max tanh(a/i5 max ) which preserves 
the sign of a but saturates at constant value, E max . 

Summing up, our algorithm requires calculating g and 
/ at each time step and choosing a field according to 



E(t) = E max t<mh(S(t)(f ■ g)/E 

max / ' 



(10) 



In order to begin the control process it is necessary to 
seed a small fraction of population into the excited states 
at initial time. This can be done by exciting the system 
with a weak pulse tuned to the optical transition. Alter- 
natively, one can begin with an initial condition slightly 
rotated from the identity. As the seeding is negligibly 
small, the emerging field will produce the target unitary 
transformation also when applied to the 'pure' initial con- 
dition U(0) = I as required. 
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FIG. 2: Schematic picture of Na2 molecular model consisting 
of ground (X E^) and excited electronic potential 

energy curves. A submanifold of the ground state levels con- 
stitutes the quantum register. The remaining levels serve as 
mediators which couple the register levels through interaction 
with the electric field. The goal is to produce a unitary trans- 
formation on the register levels with no intermediate popula- 
tion of the mediating states. 



We apply the local control technique described above 
to the implementation of an SU(N r ) Fourier transform 
on the ground potential surface of a two-electronic sur- 
face model of an Na2 molecule. The first N r vibronic 
states of the ground electronic surface, X 4 Y,+ , comprise 
the quantum register. All ground states are coupled via 
the dipole coupling to the vibronic levels of the electroni- 
cally excited surface, A 1 ^^ (see figure|5J). The final time 
for the implementation was taken to be 128 picoseconds. 
Figure |3 summarizes our results for producing an SU(6) 



Fourier transform, 
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= ^T = e 2 "/ 6 , 



on a submanifold of the two potential surfaces containing 
seven ground state vibrational levels and three excited 
state vibrational levels. The top plot shows the elec- 
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FIG. 3: Top: The control field obtained by the algorithm. 
Middle: Controlled evolution of the elements of the unitary 
propagator (absolute values). Bottom: Evolution of objective 
\fj (solid) and constraint C (dashed) showing monotonic in- 
crease in the former while keeping the latter constant. 

trie field obtained by the local control procedure. The 
middle plot shows the evolution of the absolute value of 
the unitary propagator elements under the derived field. 
The monotonic increase of the objective \fj towards its 
maximal value six can be seen in the bottom plot. The fi- 
delity of the gate, F = J/N 2 , achieved at the final time, 
is very close to unity such that log 10 (l — F) —1.7. 
Also note that only negligible leakage out of the quan- 
tum register has occurred throughout the process, which 
is apparent in the flat constant value of the constraint 
C along the evolution. A more complete picture of the 
unitary propagator elements can be obtained from fig- 
ure 0] The elements are shown in plots a), b) and c) at 
times t — 0, t ~ T/2 and t = T respectively as vectors 
in the complex plane. At time t = there are six vec- 
tors pointing along the real axis towards unity. These 
are the diagonal elements of the identity. The remaining 
elements are zero. The slight rotation of the diagonal el- 
ements and small population of the offdiagonal elements 
at initial time are due to the small seed rotation induced 
on the initial condition to initiate the control algorithm. 
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FIG. 4: The elements of the unitary propagator at times a) t = 0, b) t = T/2 and c) t = T as vectors in the complex plane. 



We stress, however, that due to the smallness of the per- 
turbation, once the field is obtained it can be applied to 
the desired 'pure' initial condition U(0) = I, with neg- 
ligeble deviations from the current results. As time pro- 
gresses the elements rotate and shrink/expand to obtain 
the SU (6) Fourier transform at the final time. It can be 
clearly seen that the elements at final time (c) divide the 
circle into six equal parts and thus are, up to an unim- 
portant phase rotation, just powers of w, the sixth root 
of unity appearing in FT(6). 

Some understanding of the mechanism by which the 
leakage is avoided can be obtained by looking at the spec- 
trum of the field E(t) and its square |-E(i)| 2 (top and 
bottom of figure [S] respectively) . The spectrum E(lu) of 




Frequency [a.u 



FIG. 5: The spectrum E(td) of the field with one-photon 
transition frequencies marked by (red) circles (top) and the 
spectrum of the field intensity |_E| 2 (o;) with two-photon tran- 
sition frequencies marked by (green) triangles (bottom). 

the field, corresponds to one-photon processes. A glance 
at the plot of E(ui) (top of figure [5j reveals that there 
are 'holes' in this spectum at precisely the points corre- 



sponding to transitions from the quantum register to the 
excited states, as indicated by the circles. These 'holes' 
are evidently responsible for the absence of excitations. 
As one-photon interactions are suppressed it must be two 
(and higher) photon processes which produce the evolu- 
tion towards the desired target. The spectrum |_B| 2 (w) 
of the field intensity corresponds to two-photon processes 
produced by absorption and immediate emission of a pho- 
ton. The bottom plot of figure [5] shows that |I£| 2 (cj) dis- 
plays peaks at the precise frequencies corresponding to 
energy differences between the register states, indicated 
in the figure by green circles. This implies that the field 
is in two-photon resonance with transitions correspond- 
ing to the register states but is detuned from one-photon 
resonance with excited state transitions. 

The model system used above illustrates the general 
features of our local control method; however it suffers 
from the following drawbacks. First, the analysis de- 
scribed above was performed in the interaction picture 
where the Hq motion is transformed away. In practice, 
however, this drift motion must be taken into account. 
Second, it is generally expected that the computational 
power scale exponentially with the physical resources 
such that for example N r physical bits carry 2 Nr entities 
of information. In the molecular system studied here, 
as in other studies of molecular quantum computation, 
this requirement is not satisfied since N r physical levels 
correspond to only N r information entities (log 2 N r bits) 
namely the scaling is only linear. However, the control 
method we have proposed does not depend on the specific 
model studied; therefore it can be applied to alternative 
models which are both driftless and scalable. 

In summary, we have shown that it is possible to sys- 
tematically design fields to produce arbitrary unitary 
transformations while avoiding leakage from the quan- 
tum register. 
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